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We present numerical simulations of coupled Ginzburg- 
Landau equations describing parametrically excited waves 
which reveal persistent dynamics due to the occurrence of 
phase slips in sequential pairs, with the second phase slip 
quickly following and negating the first. Of particular interest 
are solutions where these double phase slips occur irregularly 
in space and time within a spatially localized region. An ef- 
fective phase diffusion equation utilizing the long term phase 
conservation of the solution explains the localization of this 
new form of amplitude chaos. 
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Spatio-temporal chaos poses one of the important chal- 
lenges in the investigation of spatially extended dynam- 
ical systems. A central question regarding these sys- 
tems concerns the characterization and classification of 
the chaotic states. This has been studied most exten- 
sively for chaotic traveling waves in one dimension within 
the framework of a complex Ginzburg Landau equation. 
There, two different dynamic regimes have been iden- 
tified, a phase-chaotic regime and an amplitude-chaotic 
regime The latter is characterized by the occur- 

rence of phase slips during which the amplitude of the 
wave goes to zero and the total phase of the system 
changes by 2ir, i.e. a wavelength is inserted or elimi- 
nated. In the phase-chaotic regime essentially no phase 
slips occur |||] . This chaotic regime can be described by 
an equation for the phase of the wave alone j5| . In the 
amplitude-chaotic regime this is not possible since the 
phase equation breaks down during phase slips. 

In this letter we present a different class of spatio- 
temporal chaos, which we have identified in one- 
dimensional simulations of parametrically excited waves. 
The observed chaos is characterized by the occurrence 
of double phase slips which effectively preserve the to- 
tal phase. Thus, although the phase description breaks 
down during each of the phase slips, the long-time dy- 
namics can still be described by an effective phase equa- 
tion. The second striking result of our simulations is 
the localization of spatio-temporal chaos in a part of the 
homogeneous system. That is, there are solutions that 
include a spatial region, bounded by quiescent regions on 
each side, in which the dynamics are irregular in space 
and time. Experimentally, similar phenomena have been 
observed in Taylor vortex flow || , Rayleigh-Benard con- 
vection , and parametrically excited surface waves || . 
So far, the localization mechanism in these systems is 
only poorly understood. 



For the system of equations discussed in this letter 
(which model parametrically excited waves) the localiza- 
tion can be understood using the effective phase equa- 
tion. We confirm the validity of this approach by ex- 
plicitly determining the effective phase diffusion coeffi- 
cient. The short-scale chaotic dynamics amount to fluc- 
tuations on the large scales captured by the phase de- 
cription. They may therefore introduce a noise term 
in the effective phase diffusion equation similar to that 
arising in the Kardar-Parisi-Zhang equation ||, which 
describes the large-scale behavior of spatio-temporal 
chaos in the Kuramoto-Sivashinsky equation fll0|-|l2|l , the 
leading-order phase equation for the complex Ginzburg 
Landau equation ||. 

Parametrically excited waves arise quite generally in 
systems exhibiting weakly damped modes that are oscil- 
latory in both space and time, when these modes are 
forced at twice their natural frequency |l3|Jl^ ]. Fara- 
day's experiment in which surface waves are formed on a 
liquid in a vertically oscillated container is perhaps the 
best known example (e.g. fI1||l6| f| ) . Other examples are 
parametrically forced spin waves |l^ff^], and traveling 
waves in electroconvection jl9| and in Taylor Dean flow 
pO| . For small forcing amplitude, these systems can be 
modeled by a pair of coupled Ginzburg-Landau equations 

d T A + sd x A = dd\A + aA + bB 

+cA(\A\ 2 + \B\ 2 )+gA\B\ 2 (1) 
d T B - sd x B = d*d%B + a*B + bA 

+c*B(\A\ 2 + \B\ 2 )+g*B\A\ 2 . (2) 

The dependent variables A and B are complex, and rep- 
resent amplitudes of left and right traveling waves which 
are summed together to yield the solution of the under- 
lying system: 



u(x,t) = eA{X,T)e i{q " x - 
+c.c. + 0(e 2 ). 



-0 



eB(X,T)e 



i(q c x + ^t) 



(3) 



They vary on slow time and space scales, T = e 2 t and 
X = ex respectively. The coefficient b of the linear cou- 
pling term gives the amplitude of the periodic forcing. 

Stationary solutions to equations (1,2) of the form 
A = A e iqX , B = B e lqX (where A and B are com- 
plex constants) include those where the amplitudes of A 
and B are equal, corresponding to standing waves in the 
underlying system which are phase-locked to the exter- 
nal forcing |13|,0. In this letter we discuss results of 
numerical simulations of the ensuing dynamics when the 
standing waves are perturbed. These simulations were 
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performed using a Fourier spectral method in space (pe- 
riodic boundary conditions) and a fourth order Runge- 
Kutta/integrating- factor scheme in time. 

The linear stability diagram for the standing-wave so- 
lution of equations (1,2) with the parameters chosen for 
this study is shown in figure 1 For small values of 

b (b = 0.1, say) solutions to equations (1,2) behave like 
those of a single real Ginzburg-Landau equation; given 
an initial condition slightly to the right of the Eckhaus 
stable region, the solution will undergo a phase slip which 
reduces the wave number and moves the solution into the 
stable region (e.g. |22|]). This single-phase-slip behavior 
is shown in figure 2a. For larger values of b however 
(b = 0.6, say), the behavior of the solution is markedly 
different; given an initial condition slightly to the right 
of the Eckhaus stable region, the solution undergoes a 
phase slip which reduces the average wave number but a 
short time later undergoes a second phase slip at essen- 
tially the same location which restores the wavenumbcr 
to its original value. This double phase slip, shown in 
figure 2b, causes the solution to remain in the Eckhaus 
unstable region, thus allowing persistent dynamics. 
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FIG. 1. Stability diagram for parametrically driven stand- 
ing waves with a = —0.05, c = —1 +4i, d = 1 + 0.5i, s = 0.2, 
g = -1- 12i. 

Simulations of a small system (5 wavelengths long) 
show that the complexity of the ensuing dynamics is de- 
pendent on the distance of the initial condition from the 
Eckhaus stable region. Very close to the Eckhaus curve 
simple periodic solutions occur. In these solutions dou- 
ble phase slips occur at only one location in space, pe- 
riodically in time. Further to the right, away from the 
Eckhaus stable regime the simple periodic solutions lose 
stability to more complicated periodic solutions and these 
more complicated periodic solutions in turn lose stability 
to solutions where the double phase slips occur irregu- 
larly in time and space over the entire domain. This 
irregular occurrence of double phase slips is a form of 
spatio-temporal amplitude chaos. Starting with an ini- 
tial condition far to the right of the Eckhaus curve, single 



phase slips occur which reduce the wavenumber and thus 
move the system to the left, into a region where a sta- 
tionary, periodic, or chaotic solution is stable. 
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FIG. 2. Space-time diagrams illustrating a) single phase 
slip observed for small values of the forcing amplitude (e.g. 
b — 0.1), and b) double phase slip observed for larger values 
of the forcing amplitude (e.g. b = 0.6). Other parameters as 
in figure 1. 

In larger systems (here 47 wavelengths long) and at in- 
termediate distances from the Eckhaus stable regime we 
find that the chaotic activity need not spread through 
the whole system. Instead, double phase slips occur only 
in a spatially confined region. This is illustrated in figure 
3 which shows the location in space and time of dou- 
ble phase slips. Here an initial maximum in the local 
wavenumber triggers a double phase slip which subse- 
quently causes more double phase slips in the neighbor- 
hood of the first. These double phase slips occur irregu- 
larly in space and time within a region whose size initially 
grows but remains bounded in space. Outside of this 
region the solution remains stationary (i.e. in physical 
space there are regular, periodic standing waves, cf.(3)). 
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FIG. 3. Space-time diagram showing each double phase slip 
as a single dot. The spatio-temporal chaos is confined to a 
small part of the system (6 = 0.6, other parameters as in figure 
1; the averaged effective wavenumber is < q >= 0.377). 

To understand the localization mechanism we turn to 
the phase diffusion equation 



d T cj>{X,T) = D{q)d x 4>{X,T), 



(4) 



which describes the slow evolution of the phase (j> of a 
steady pattern. Equation (4) is nonlinear through the 
dependence of the diffusion coefficient on the wave num- 
ber q = d x <p. The use of the phase equation is motivated 
by the observation that the temporally averaged local 
wavenumber is larger in the chaotic region than in the 
steady region. The state shown in figure 3 can therefore 
be considered to consist of two domains, one with large 
and one _with_small wavenumber. It has been shown pre- 
viously 
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|2J] that such domain structures can be sta- 
ble if the diffusion coefficient D(q) is negative (indicating 
an unstable solution) only over a small range of wavenum- 
bers q. Initial conditions with a uniform wavenumber in 
this range will then evolve to a structure consisting of 
domains in which the local wavenumber lies in either of 
the two adjacent ranges of stable wavenumbers. The sta- 
bility of the domain structure is due to the conservation 
of the total phase and the instability of the state with 
the corresponding uniform wavenumber. 

As a basic condition for this mechanism to explain the 
localization shown in figure 3 the phase description has 
to be valid, requiring in particular that the phase is con- 
served. Across individual phase slips the phase is not 
conserved and the phase equation (4) breaks down. How- 
ever, across a double phase slip the phase is conserved. 
On sufficiently long time scales we therefore expect an 
effective phase equation to be appropriate. Reflexion 
symmetry in space but not in time suggests an equation 
of the form (Q) with the phase and the diffusion coeffi- 
cient replaced by an effective phase <fi an d an effective 
diffusion coefficient D(q). To substantiate this claim we 
determine D(q) by measuring the (diffusive) response of 
the extended chaotic state to a localized time-periodic 
forcing. Specifically, we include a spatially local, time- 
periodic advection term in equations (1,2) giving: 



d T A + (v + s)d x A = dd x A + aA + bB 

+cA(\A\ 2 + \B\ 2 )+gA\B\ 2 
d T B + (v- s)d x B = d*d 2 x B + a*B + bA 



(5) 



with v 



c*B{\A\ 2 + \B\ 2 ) +g*B\A\ 2 (6) 

v sin(cjT) for Xy < X < X 2 
otherwise. 



For relatively small values of wo formerly stationary so- 
lutions start to drift within the region X\ < X < X 2 , pe- 
riodically reversing their direction according to the sign 
of v. For the regions X < X\, and X > X 2 the situation 
resembles an imposed boundary condition at X = X%, 



and X — X 2 at which the phase of the solution varies 
sinusoidally. This is similar to the approach used in ex- 
periments on turbulent Taylor vortex flow wherein an 
endcap was moved sinusoidally fl25|. 
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FIG. 4. Decay of the magnitude of the coefficient of the 
time Fourier transform of A corresponding to the frequency 
ljo in the localized advection term v in equations (5,6), for two 
different values of ujq (b — 0.61, other parameters as in figure 
1; the averaged effective wavenumber is < q >= 0.4006). 

The solution to the diffusion equation dt<j>(X,t) — 
Dd%4>(X,t) for X > X 2 , with 4>{X 2 ,t) = 4> sm(Lu t), 



4>(X, t) = 4> e- a ( x - x ^ sm(LU t - (3(X - X 2 )) (7) 
2D 



where a — [3 



(8) 



If the phase of the solution A in the region X > X 2 obeys 
this relationship, then the magnitude of the Fourier mode 
A(X, wo) corresponding to the frequency luq will decay 
exponentially in space as X increases from X 2 . Figure 
4 shows that this exponential behavior is realized in our 
simulations. Note that the data for figure 4 is from a 
simulation that produced double phase slips irregularly 
in space and time over the entire domain (not localized 
chaos). The phase slips introduce noise which can be 
seen in figure 4. In order to reduce the effect of this noise 
and extract a reliable decay rate a, the Fourier integral 
has been extended over many periods (88) of the function 
v{X, T) in the advection term in equations (5,6). Equa- 
tion (8) shows that for diffusive behavior the decay rate 

a is proportional to ^J 2 . The ratio of the slopes of the 
curves in figure 4 are consistent with this diffusive scal- 
ing. We conclude that the dynamics are indeed diffusive 
and show in figure 5 the diffusion coefficient as a function 
of the wavenumber. The solid curve shows the analyti- 
cal phase diffusion coefficient for the stationary solution 
pl| while the triangles show the effective phase diffusion 
coefficient for the chaotic solution as measured using the 
decay rate a of the Fourier mode A(X, uj) for ui = ujq and 
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lu = —luq for X > X2 and X < X\ as described above. 
Despite the scatter in the data it is clear that D(q) de- 
creases with decreasing q and presumably goes to zero at 
a wavenumber for which the diffusion coefficient of the 
nonchaotic state is still negative as shown by the dashed 
line. Thus the system is diffusively unstable over a range 
of wavenumbers, and initial conditions with an average 
wavenumber in this range will evolve into domains of 
chaotic and nonchaotic waves. 
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FIG. 5. Analytical phase diffusion coefficient for the sta- 
tionary solution (solid line) and effective diffusion coefficient 
for the chaotic solution (triangles represent numerical results, 
the dashed line is the presumed extension of these results). 
(b = 0.61, other parameters as in figure 1.) 

Figure 5 suggests that the chaotic activity should not 
disappear homogeneously if the average wavenumber is 
decreased towards the stable regime. Instead, the homo- 
geneously chaotic state should split up into chaotic and 
stationary domains in which the local wavenumber is in 
the respective stable regimes (D > and D > 0). This 
is indeed found in our numerical simulations. 

To summarize, we have presented numerical simula- 
tions of spatio-temporal chaos in parametrically excited 
standing waves. We have demonstrated that despite the 
fact that the chaos must be classified as amplitude chaos 
its large-scale features can nevertheless be described by 
an effective phase equation. The reason for this surpris- 
ing result lies in the fact that the phase slips always come 
in pairs such that the total phase is effectively conserved. 
We have determined the effective diffusion coefficient. 
This allowed us to identify a simple mechanism for the 
localization of the spatio-temporally chaotic state which 
is related to the occurence of domain structures and zig- 
zags in steady patterns . 

So far we have not identified a particular physical sys- 
tem which exhibits the behavior described here. It is 
natural to suspect that the closing of the Eckhaus curve 
(cf. figure 1) is related to the occurence of double phase 
slips. If that should be the case, parametric driving of 
waves which are Benjamin- Feir unstable in the absence of 



any forcing should be a good candidate for this behavior. 
Within the framework of the coupled Ginzburg-Landau 
equations (1,2) these waves become unstable near the 
band center when forced sufficiently strongly j27j. De- 
pending on the sign of certain nonlinear coefficients this 
indicates instability at all wavenumbers and a closing of 
the Eckhaus curve. 
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